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ABSTRACT 

We present a spectroscopic study of the most massive stars in the young (4 Myr old) 
stellar cluster LH 95 in the Large Magellanic Cloud. This analysis allows us to complete 
the census of the stellar population of the system, previously investigated by us down 
to 0.4 solar masses with deep HST Advanced Camera for Surveys photometry. We 
perform spectral classification of the five stars in our sample, based on high resolution 
optical spectroscopy obtained with 2.2m MPG/ESO FEROS. We use complementary 
ground-based photometry, previously performed by us, to place these stars in the 
Hertzsprung-Russel diagram. We derive their masses and ages by interpolation from 
evolutionary models. The average ages and age spread of the most massive stars are 
found to be in general comparable with those previously derived for the cluster from 
its low mass PMS stars. We use the masses of the 5 sample stars to extend to the 
high-mass end the stellar initial mass function of LH 95 previously established by 
us. We find that the initial mass function follows a Salpeter relation down to the 
intermediate-mass regime at 2 M . The second most massive star in LH 95 shows 
broad Balmer line emission and infrared excess, which are compatible with a classical 
Be star. The existence of such a star in the system adds a constrain to the age of the 
cluster, which is well covered by our age and age spread determinations. The most 
massive star, a 60-70 M Q 02 giant is found to be younger (< 1 Myr) than the rest 
of the population. Its mass in relation to the total mass of the system does not follow 
the empirical relation of the maximum stellar mass versus the hosting cluster mass, 
making LH 95 an exception to the average trend. 

Key words: stars: early-type; stars: massive; stars: luminosity function, mass func- 
tion; Magellanic Clouds; open clusters and associations: individual: LH 95 



1 INTRODUCTION 



It is well established I Lada fc Lada|2003 1 that the vast ma- 
jority of stars are not formed in isolation, but as part of 
stellar associations and clusters. These comprise hundreds 
or even thousands of stars along the entire mass spectrum, 
from the most massive OB-type down to low-mass late-type 
stars. The study of young (few Myr old) clusters, there- 
fore, provides a fundamental tool to understand how the 
star formation process takes place. Both high- and low-mass 
stars are, from different points of view, of complementary 
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importance in such clusters. On one hand, the short-lived 
massive stars are the direct evidence of the youth of these 
systems, and their energy output through winds and UV 
radiation is responsible for the gas removal, a process that 
strongly affects the survival of the stellar system. On the 
other hand, the low-mass stars constitute the majority of 
the newly formed stellar population, and their slow evolu- 
tion during the pre-main sequence (PMS) phase traces the 
timescale of star formation. 

The study of star-forming regions in the Milky Way 
is affected by the high extinction produced by the dust 
in the Galactic disk, limiting the depth of optical obser- 
vations. This problem does not affect significantly studies 
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in the Magellanic Clouds (MCs), our nearby metal-poor 
dwarf galaxies. Both the Large and Small Magellanic Clouds 
(LMC, SMC) are characterized by an exceptional sample of 
star-forming regions, the Hn regions (|Henize||1956 Davies 



L 

19761, where hydrogen is already being ionized by the 



UV-winds of massive stars. A plethora of young stellar as- 
sociations and clusters are embedded in these regions (e.g., 



iBica & Schmitt 1995; Bica et al. 1999), giving evidence of 



current clustered star formation 

LH 95 (|Lucke k. Hodge|1970|) is a young stellar cluster 



in the LMC related to the Hn region LHA 120-N 64 (Henize 



1956), located to the north-east of the super-bubble LMC 4. 
The association has a projected size of ~ 10 x 20 pc, and is 
not centrally concentrated but presents at least 3 main sub- 
clusters. Its stellar population has been studied by |Goulier^l 
mis et al. ( 2002 1 using ground-based BVR and Ha photom- 



etry; they determined the initial mass function (IMF) from 
30 Mq to 2 M©. More recently, very deep Hubble Space 
Telescope Advanced Camera for Surveys (HST/ACS) ob- 



servations in _F555W and F8UW bands (Gouliermis et al 



2007[ ) led to the disc overy of a populati on of ~ 2500 PMS 
stars in the region. In Da Rio et al. ( 2009 hereafter, Paper I), 
using these data, we derived the first extragalactic IMF in 
the subsolar regime, complete down to ~ 0.4Mq while mem- 
bers have been detected with masses as low as O.2M0. Later, 
Da Rio et al. (2010[ hereafter Paper II), we investigated 



the age and age spread of the intermediate- and low-mass 
stars in LH 95 through a statistical technique that accounts 
for physical and observational biases in the observed colors 
and magnitudes. We found an average age of the system of 
the order of 4 Myr, and an age spread well represented by a 
gaussian distribution with a standard deviation of 1.2 Myr. 

Despite the large luminosity (and mass) range spanned 
by our ACS observations, our stellar sample was incomplete 
to the high-mass end, since five of the brightest stars were 
saturated in the HST images. Moreover, optical photometry 
of high-mass stars (with M ^ 8 M©) alone is not sufficient 
for the accurate derivation of their stellar parameters. This 
is because these stars, due to their high effective temper- 
atures T c ff, emit a large fraction of their flux in the UV 
wavelength range. As a consequence, the optical bands are 
in the Rayleigh regime and the observed colors tend to be 
independent of T e g . This well-known problem was made ev- 
ident also in Paper II, where we demonstrated that dating 
massive stars based only on their photometry is highly un- 
certain. 

Stellar parameters of massive stars can be determined 
accurately through spectroscopy. Therefore, in order to 
study the brightest stellar members of LH 95, and to con- 
struct the H-R diagram of its massive stellar population, we 
perform here a spectroscopic follow-up to our previous stud- 
ies. We acquired optical spectra of the five brightest stars 
in the region and we complemented these data with exist- 
ing BVR photometry of LH 95, obtained from the ground 



(Gouliermis et al. 20021. The aim of this investigation is 



twofold: 1) the completion of the IMF derived in Paper I 
toward the highest masses, and 2) the analysis of the ages 
of the high-mass members (and possible variations), in rela- 
tion to the similar results we presented in Paper II for the 
low-mass population. 

This paper is organized as follows. The spectroscopic 
observations and their data reduction are described in Sec- 



tions [2] and [3] In Section H] we present the classification of 
the sources. In Section [5] we derive the Hertzsprung-Russel 
diagram for the high-mass end of the LH 95 population, and 
assign ages and masses to the sources. A discussion on these 
results is presented in SectionlSJ and in SectionlTlwe discuss 
the high-mass end of the IMF. We then focus in Section 
[8] on the brightest source in the sample, where we present 
evidence that classifies it as a classical Be star. Finally, in 
Section [9] we summarize the main findings of this paper. 



2 THE OBSERVATIONS 

This spectroscopic study focuses on the brightest stellar 
sources in LH 95, which were saturated in our HST/ACS 
images. For the photometry of these sources we use the 



measurements by Gouliermis et al. (20021, based on BVI 



ground-based observations. While this photometric study 
achieved significantly shallower detection limits and poorer 
spatial resolution than our ACS imaging, it provides accu- 
rate BVI photometry for all the bright sources within the 
star-forming region. Based on this study we excluded from 
our analysis the bright stars with significantly red colors, 
as probably being field RGB stars. This is because at the 
age of LH 95 (~ 4 Myr), the turn-off from the main se- 
quence occurs at a large mass (~ 40 M©), and the RGB 
stars above this mass are not expected to have tempera- 



tures lower than 6000 K (Girardi et al. 2002 Schaerer et 



al.|1993l, thus should not show very red optical colors. The 



remaining five bright main-sequence stars are the targets of 
this spectroscopic follow-up. We annotate these sources with 
the letters A, B, C, D, and E in the color-composite VI im- 
age of our HST/ACS observations, shown in Figure n] (left 
panel). The B — V,V color-magnitude diagram (CMD) from 
the photometric study of |Gouliermis et al. (20021 is shown 
also in Figure [j] ( right panel). In this plot stars identified by 
Gouliermis et al. ( 2002 1 in a significantly larger field-of-view 



(FOV) than that of ACS are plotted with grey symbols. The 
stars observed within our ACS FOV ( [Gouliermis et al 2007) 
are highlighted with large black symbols, and the targets of 
this spectroscopic study with red symbols. They are also 
annotated in this CMD with their alphabetic IDs. 

The observations have been carried out with the 
Fibre-fed, Extended Range Echelle Spectrograph (FEROS), 
mounted on the 2.2 m MPG/ESO telescope located at the 
La Silla Observatory in Chile. This instrument has a resolu- 
tio n of R = 48, 000 and a wa velength coverage of 3600 - 9200 
A (Kaufer & Pasquini 1998 1. The spectrograph is fed by two 



fibres with an aperture of 2.0 arcsec, providing simultane- 
ous spectra of a stellar target, plus either the sky or one of 
two calibration lamps (wavelength calibration and flat-field). 
Since for our purposes the wavelength calibration is not crit- 
ical, we used the "OBJSKY" instrument setup, acquiring a 
position in the sky in parallel to the targets separated by 
2.9 arcmin. The necessary calibration acquisitions have been 
carried out separately for each observing night. Observations 
were carried out in two epochs: November 2008 and January 
2009. The summary of the observations is given in Table [l] 
Multiple spectra have been taken for the same targets to en- 
able cosmic ray detection and removal. The overall exposure 
times have been higher for fainter sources in order to reach a 
sufficient signal-to- noise ratio (SNR). Specifically, our final 
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Figure 1. Le/t panel: The five targets selected for our spectroscopic follow-up, onto the color-composite image of LH 95 from ACS 
observations (Image credit: NASA, ESA and D. A. Gouliermis - MPIA). Right panel: CMD from the |Gouliermis et al.| ( |2002^ photometry. 
The black dots indicate sources located within the ACS FOV. The five targets selected for spectroscopy are marked with thick red dots. 



Table 1. Summary of the observations 



ID ID2 



(GO20 



alternative nam 



RA Dec V exposure time 

(J2000.0) (J2000.0) [mag] [s] 



airmass 



epoch 



13 


9141 


C 


9143 


D 


18 



9136 SK -66 170 05:36:58.945 -66:21:16.130 13.17 

SK -66 172 05:37:05.553 -66:21:34.950 13.58 

SK -66 174 05:37:15.723 -66:21:38.355 13.74 

05:36:58.007 -66:21:42.613 14.53 

85 05:37:15.129 -66:21:44.304 15.58 



3600 


1.31 


Nov 2008 


1000 


1.35 


Nov 2008 


1500 


1.39 


Nov 2008 


1200 


1.26 


Jan 2009 


1200 


1.25 


Jan 2009 


1200 


1.30 


Nov 2008 


900 


1.32 


Nov 2008 


900 


1.27 


Nov 2008 


900 


1.28 


Nov 2008 


1000 


1.25 


Nov 2008 


1000 


1.26 


Nov 2008 


1800 


1.26 


Jan 2009 


1800 


1.28 


Jan 2009 


1400 


1.28 


Nov 2008 


1400 


1.31 


Nov 2008 


1400 


1.34 


Nov 2008 


1800 


1.31 


Jan 2009 


1800 


1.35 


Jan 2009 



a catalog numbers of Gouliermis et al. (2002) 
Sanduleak (1970) identification numbers 

spectra reach an average SNR in the optical range spanning 
from ~ 100 for star A, to ~ 30 for star E. 



3 DATA REDUCTION 

The basic data reduction of the spectra was performed by 
using the automatic online pipeline. This pipeline extracts 
the 39 cross-dispersed spectral orders and combines them to- 
gether in a single, wavelength calibrated, spectrum. We man- 
ually subtracted the relevant sky emission from each spec- 
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Figure 2. Reduced spectra for the five targets, covering the violet-blue range of the spectrum, relevant for spectral classification of early- 
type stars. The main stellar lines we have identified in the spectra are indicated as well. The features identified in the spectra are: He 1 
AA4009, 4026, 4121, 4144, 4387, 4471, 471 3; He II AA4200, 4541, 4686; C III/O II A4070; Si IV AA4089-4116; Si III AA4552- 4568-4575; 
N II A3995; N III AA4097, 4634 - 4640 - 4642; N IV A4058; N V AA4604 - 4620; O II A4415; C III AA4647 - 4650 - 4652. 



trum; then we normalized the continuum of each spectrum 
to unity by fitting a high-order polynomial to the continuum 
smoothed with a median filter 12A in size. Each derived con- 
tinuum has been visually inspected together with its origi- 
nal spectrum to ensure a correct estimation. The continuum 
normalization also allows us to correct for shape variations 
in the spectra of the same star introduced by changes in 
airmass and atmospheric extinction, as well as flux factors 
due to different exposure times. Cosmic rays are detected 
as > 3cr peaks above the continuum on each spectrum and 
absent in the other spectra of the same star. The average 
local noise, a, was computed in a 2 A neighborhood of each 
wavelength point. The cosmic rays were then masked, and 
the spectra were combined in one final spectrum with the 



use of a weighted average method, accounting for the in- 
dividual exposure times. In order to facilitate the spectral 
classification, the original spectra were also degraded to a 
resolution of 1A. 



The blue wavelength range of the final spectra are pre- 
sented in Figure [2] In this plot the spectra are blue shifted 
by an amount 8X/X — 10 -3 , to correct for the radial velocity 
of the system, v r — 302.8 km s~ , determined from the lines 
in the spectra. For the four faintest stars, we have manually 
removed part of the spectra at 4695A ^ A ^ 4720A, where 
an artificial feature in emission dominates the continuum of 
the 16th spectral order. By visual inspection on the individ- 
ual spectral orders, we verified that the He II A4686 line is 
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not affected by this artifact, being located at the red end of 
the 15th order. 



4 SPECTRAL CLASSIFICATION 

We assign spectral types to our targets by comparing their 
spectra to stellar atlases available in the literature. In par- 
ticular we consider the atlases of |Walb orn & Fitzpatrick 
(19901 and Sota et al. (20111, which contain a large selec- 



tion of early type stars and provides a fine coverage of the 
spectral types and luminosity classes of OB stars. We also 



use the atlases of OB type stars of Walborn et al. ( 1995 1 



and |Walborn et ah] ( |2002[ ). Based on our classification, we 
assume an uncertainty in the assigned spectral types of half 
a spectral subtype for stars A to D and one subtype for star 
E. 



Star A 

This is the brightest stellar source in the region (V ~ 
13.5 mag; Gouliermis et al.|2002 l. Due to the long total ex- 
posure time allocated to this source (exceeding 2 hours) the 
derived spectrum reaches a very high SNR (> 100). The line 
ratio between He I A4471 and He II A4541, together with the 
ratio Si III A4552/Si IV A4089 indicate a BO. 2 type star. As 
far as its luminosity class is concerned the spectrum of this 
star shows some inconsistencies. In particular, while the sig- 
nificant strength of the Si IV lines suggests a very luminous 
star of class I or II (also considering the low metallicity of 
the LMC), the presence of the He II A4686 line is consistent 
only with a luminosity class IV or V. Therefore, we assign 
to this source a giant luminosity class and we classify this 
star as a BO. 2 Hip. As we discuss later (see Section Is!) this 
star shows prominent Ha and H/3 emission, behaving as a 
classical Be star. 



Star B 



This star was already classified by Walborn et al. ( 1995 1 as 
an 03 lll(f*), and then reclassified in|Walborn et al.| (|2002 1 



as 02 III(f*), with evidence in support of a binary sys- 
tem with a later OB-type star. Our spectrum confirms this 
classification, mainly based on the prominent N IV A4058 
emission, the NV AA4604 — 4620 absorption, and the ab- 
sence of N III A4640 emission. The significant strength of the 
N IV A4058 line indicates a luminous star (^ III), but the 
He II A4686 in absorption excludes a supergiant. Accord- 
ing to Walborn et al. (20021, the spectrum is a composite 



including an unresolved later-type star, as revealed by the 
detection of the He I A4471. In our spectrum (see Figure p| 
we also clearly detect He I A4387, which did not appear in 



Walborn et al. ( 1995 1 spectrum. This further supports the 



claim that this source is a composite system. For an addi- 
tional discussion on the composite nature of this source, we 
refer to Section fe.l.ll 



Star C 

According to the ratio between He II A4541 and He I A4471 
(which are equally strong for 07 stars) this source is an 07.5 
type. The presence of N III AA4634 - 4640 - 4642 emission 



is typical of a luminous star of this spectral type, and their 
strength suggests at least a 11(f) luminosity class. The strong 
N III A4097 absorption is also supportive of a bright giant. 
We also notice a modest emission in the C III AA4647 — 
4650 — 4652, which is atypical for this spectral type. The 
strength of the He II A4686 absorption in contrast with the 
N Illemission, indicating a luminosity not higher than class 
III. We argue that the unusual strength of the N emission 
lines might be due to nitrogen enrichment, or, alternatively, 
could be indicative of a composite nature for this source. As 
a compromise, we classify this star as 07.5 111(f). 

Star D 

The ratio between He II A4541 and He I A4471 indicates a 
type of 06.5 for this star. In addition, the presence of strong 
He II A4686 absorption, typical of main sequence O-type 
stars, leads to the classification of this star as an 06.5 V. 



Star E 

Given the strength of the He I and the absence of He II lines, 
this source is an early B-type star. The strength of Si IV and 
Si III, together with the absence of Si II AA4128 - 4130 
indicate a B1.5 star. The additional presence of the N II 
AA3995 suggests a luminous star, but not a supergiant given 
the absence of several indicators of such high luminosity 
(e.g., strong O II lines). Therefore we assign the spectral 
type B1.5 III. 



5 THE H-R DIAGRAM 

In this section we derive the stellar parameters of the ob- 
served sources from both our spectroscopic analysis de- 
scribed above, and the available photometric measurements 
of these sources ( |Gouliermis et aLJ|2002| |, and we construct 
the H-R diagram for these stars. We consider appropriate 
temperature scales for early-type stars with LMC metallic- 
ity ([Fe/H]~ -0.3; |Luck et aL]|1998| >, and we derive the 
effective temperatures (T c s) of the stars from their spectral 
types. There is a number of spectral type - T c s conversions 
available in the literature, covering different ranges in the 
parameter space. The derivation of such relations usually 
relies on model stellar atmospheres, in which T e g is an in- 
put parameter. The associated spectral types are obtained 
from the presence of diagnostic lines, following observational 
criteria such as those we applied in the previous section. 

The stellar atmosphere modeling in particular for O- 
type stars is quite challenging (see, e.g., Kudritzki & Puis 



2000 1 , since it requires several factors to be considered for 



the effects of metals in the atmospheric structure and emer- 
gent spectrum, such as a non-local thermodynamic equilib- 
rium (NLTE) treatment, spherical expansion due to stellar 
winds, and line blanketing ( |Martins et al.|2005[ ). Such stud- 
ies show a clear dependence of the temperature scale on 
the luminosity class, with luminous stars being colder than 
dwarfs at a given spectral type. In addition, |Massey et al.| 
(2005) found evidence for a systematic change in the temper- 



ature scale with metallicity. They collected high SNR optical 
and UV spectra of a sample of O-type stars in the Galaxy 
and in the MCs, and derived the stellar T a a via modeling. 
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These authors found that SMC O-type stars are systemat- 
ically hotter than Galactic stars at the same spectral type, 
while LMC stars follow an intermediate temperature scale 
between the two, as expected by their metallicity. In agree- 
ment with previous works, they also found lower T e ff for 
supergiants (luminosity class I), while no significant differ- 
ences were measured between giants and dwarfs (III and V). 
Considering these results, for the late O-type stars in 
our sample we adopted the |Massey et aL| ( |2005| ) LMC scale, 
which covers the range 03 - 09.5. Temperatures are not 
explicitly provided for 02 stars, but the sample of Massey 
|et al.| ( |2005| ) includes two LMC 02 III stars, to which they 
assign a lower-limit of T c ff = 48000 K, and two 02 V with 
temperatures between 51000 K and 55000 K. Based on these 
measurements, we assign a T cff = 50000 K ± 2000 K to our 
02 III star. This value is also compatible with the Massey 
V+III scale extrapolated to 02 stars and in agreement with 
other 02 giant stars in the LMC (e.g., |Evans et al.||20lo| |. 
For the two early B-type stars, we adopt the temperature 



scale of Bessell et al. (19981, which matches very well the 



Massey scale at B0. We note that for B-type stars, the un- 
certainties in the spectra modeling become less prominent. 
This is supported by the fact that all the recent temperature 
scales for O-type stars (e.g., Vacca et al.|1996 Martins et al 



2005 Massey et al.|2005) are in agreement with each other 



at the O-B transition predicted at T c s — 30000 K, whereas 
they present significant differences at higher T e g. Moreover, 
in the B-type range we do not expect any significant depen- 
dence of the T e ff scale on metallicity, as shown by Massey et 
aLlpOOSl. 



We derive the bolometric luminosities of our stars by 
first correcting our V magnitudes for extinction, and then 
applying bolometric corrections according to the stellar tem- 
peratures. Extinctions are estimated from the E(B — V) 
color excess in the photometry of |Gouliermis et al.| (J2 002). 
We use the intrinsic (B — V)o colors published by Martins & 



Plez| (|2006) for O-type stars, and by |Bessell et al.| ( |1998| > for 
B-type stars. Since O-type stars are well into the Rayleigh- 
Jeans regime at optical wavelengths, the intrinsic colors are 



hardly sensitive on T B s ■ According to Martins & Plez ( 2006 1 , 
(B — V)o decreases from 0.28 mag to only 0.26 mag in the 
entire O spectral class, with differences between luminos- 
ity classes less than 0.01 mag. Subsequently, we also as- 
sume that metallicity variations do not affect the intrinsic 
(B — V) color. For the two early-B stars in our sample we 
consider the intrinsic colors from Bessell et al. (19981, as- 
suming log(? = 4.0. According to this study, at the Teff of 
early B-type stars, changes in surface gravity of more than 
1 dex cause very small variations of (B — V)o, of the order 
of 0.01 mag. As a consequence, the uncertainties in the lu- 
minosity classes of sources A and E are not expected to bias 
in any way the assumed intrinsic colors of the two stars, and 
thus the derived optical extinction. 

The measured color excess E(B — V) is then con- 
verted into Av assuming the average LMC reddening law 
Rv = 3.41 ± 0.06 (|Gordon et "aLl|2003[), a value slightly 



higher than the typical Galactic, Rv =3.1 (Cardelli et 
|al.||1989"| . The photometric errors in B and V are negli- 
gible, given the high luminosity of the bright members of 
LH 95. Nevertheless, we add a constant 0.035 mag uncer- 
tainty to every star in each band, representing the precision 



Therefore, the derived E(B — V) has an error of 0.05 mag. 
We convert the extinction corrected magnitudes, Vo, into 
absolute magnitudes, Mv, considering a distance modulus 
(m - M)o = 18.41 mag (Paper I). 

Finally, we convert absolute magnitudes into total stel- 
lar luminosities (log L/L@), by applying the appropriate 
V— band bolometric corrections (BCv)- For our O-type 



stars we use the relation from Massey et al. ( 2005 1 



BC = 27.99 -6.90- log T eff . 

Clearly, as for the intrinsic colors, these authors found BCs 
to be mainly dependent on temperature, and independent on 
luminosity. Also, they found good agreement with the BCs 
derived from previous works ( Vacca et al.|1996 l. A compar- 
ison of this relation with the bolometric corrections of lMar-l 
tins et al.| ( |2005[ ) yields differences below 2%. For the BCs 
of B-type stars, we consider again the values from |Bessell et| 
aLl (1998 l. 

Star A, as we describe in detail in Section [8] appears 
to be a Be star, showing prominent Ha and H/3 emission. 
This excess, however, does not affect our derived stellar pa- 
rameters for this star. The Ha line is at wavelengths longer 
than the V-band bandpass; on the contrary, the H/3 is inside 
the B— band throughput. As we show in Figure [91 however, 
this line emission appears to be modest, with an equivalent 
width (E.W.) of ~ 2. 8 A. This is small compared to the E.W. 
of the B— band filter profile, which is ~ 912A. Furthermore, 
at the H/3 wavelength (A = 486lA) the filter throughput is 
about 40% its peak. As a consequence, the line emission we 
detect for this star produces a measured flux excess in our 
B— band photometry of the order of 0.1%, or 0.001 mag. 
This is therefore absolutely negligible for the derivation of 
the Av from the E(B — V) color excess. 

We derive the bolometric luminosity from the relation 

log L/Lq = 0.4 ■ (Mboi,© -M v - BC V ), 

where Mboi,© = 4.75 is the absolute bolometric magnitude 
of the Sun. The parameters derived for our five stars in this 
section are given in columns 1 to 4 of Table [2] The positions 
of the stars in the H-R diagram, as well as the errors due to 
both T c ff uncertainties (diagonal error bars) and photometry 
(vertical error bar hats), are shown in Figure [3] In the same 
plot we overlay a set of Geneva evolutionary models for the 
LMC metallicity (Z = 0.008). These are from |Schaerer et] 



al. ( 1993 1 and are computed assuming no stellar rotation. In 



general, however, stellar rotation plays an important role in 
the evolution of massive stars ( Maeder fe Meynet|2000 |, by 
enhancing mass loss, introducing anisotropy, inducing trans- 
port of angular momentum and mixing of chemical elements 
in radiative zones. While these effects influence the stellar 



evolution in the H-R diagram (Meynet & Maeder 20001, 



they are critical mainly for the evolved stages of massive 
stellar evolution, and in particular the Wolf-Rayet phase 
(Meynet & Maeder 2003) and the final pre-supernova stages 
( [Hirschi et al.||2004| ). In Figure |3| we also plot with dashed 
lines evolutionary tracks from the grid of models by |Schaerer| 
et al. (1993[) for selected stellar masses. Unfortunately there 



of the photometric calibration of Gouliermis et al. (2002) 



are no evolutionary models considering rotation, which are 
available for the metallicity of the LMC and for the en- 
tire high-mass stellar range, except of those computed by 
|Meynet fc Maeder| ( J2005| ) for masses of 40, 60 and 120 M Q , 
which include an assumed rotation of vm = 300 km/s. We 
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Table 2. Derived stellar parameters for our target stars. 



ID 



T c ff 
(K) 



BC V 

(mag) 



(mag) 



log L/LgM masJ°l age 

(MqV (Myr) 



A 29170 

B 50000± 

C 35480 

D 38110; 

E 24380+ 



+270 
-330 

2000 
2000 

f660 
-650 

f800 
-610 

1990 

4330 



-4.4337.0 [HI 0.644 ±0.174 



5.474to.oig±o'o7 33 ' 2 ± 2 - 7 5 - 43 ± °- 41 



-,+0.047+0.07 



-3.407 



+0.060 



0.464 ±0.174 



-2.430+o4|° 0.982 ±0.174 



q+0. 023+0. 07 



5.073to.o2o±o'o7 26.8 ±1.3 2.05 ±1.14 



q+0. 136+0.07 



a The first error in log L/Lq represents the propagation of the uncertainty in log T e g (and conse- 
quently in BC), the second error represents the propagation of the uncertainty in My. 



The masses and ages of the stars are determined in Section 



5.1 



O 




4.8 



4.6 



4.4 
log T 



4.2 



4.0 



Figure 3. H-R diagram of the selected targets, together with evolutionary tracks for initial masses between 10 and 120 Mg (dashed 
lines) and isochrones for ages between 0.4 Myr and 40 Myr (solid lines) from the Geneva models of |Schaerer et ajT] ( |1993[ l for LMC 
metallicity. The dotted lines are evolutionary models for stars with rotation, for masses of 40, 60 and 120 Mq and LMC metallicity, from 
|Meynet fc Maeder| 12005) . 
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show also these tracks in Figure [3] with green dotted lines. 
For clarity we have not plotted the evolved stages for these 
masses. From the comparison between the two families of 
tracks shown in the Figure, it is evident that, at the H-R 
diagram positions of our five stars, the consideration of ro- 
tation or not does not affect significantly the derived stellar 
parameters. 

5.1 Masses and Ages of the stars 



Considering the above, we utilize the grid of Schacrer ct al. 
( 1993 1, and assign masses and ages to our stars by interpo- 



lation. A particular attention was given to the evaluation of 
the uncertainties associated to these parameters. As shown 
in Figure [3] the errors in log T e g and log L are correlated 
and asymmetric. Moreover, the probability for the star to be 
measured in different positions of the H-R diagram within 
the error may not be uniform. This is for example the case 
of source E, the measured logT c ff and logL of which are lo- 
cated over the H-exhaustion locus of the evolutionary mod- 
els, i.e., the turn-around feature of all the isochrones and 
tracks in the figure at logT ~ 4.4. After the H-exhaustion, 
the evolution of a star is much faster in the H-R diagram, 
and thus, the chance of detecting a star at colder temper- 
atures than this point is significantly lower than at hotter 
temperatures. We take into account this effect to normalize 
our results by introducing a prior to the estimation of the 
masses of our targets. 

In order to evaluate the most probable masses and ages 
for our stars, we proceed with the application of a Monte 
Carlo (MC) technique, as we describe in the following para- 
graphs. First, we construct the two-dimensional density dis- 
tribution of the relative rapidity of stellar evolution on the 
H-R diagram. For this construction we interpolate the orig- 



inal grid of isochrones and tracks of Schaerer et al. ( 1993 1 



onto a dense grid with uniform spacing in mass and age. For 
every point in the new grid, the relative speed of evolution 
is proportional to the inverse of the distance, on the same 
track, between neighboring points. This distance depends 
on the assumed metric, for which we choose a linear scale 
for both axis of the H-R diagram. Since we are interested 
in the relative variation of evolution speed between different 
points in the H-R diagram, the units chosen in both axis to 
derive the distance between points are arbitrary and equiv- 
alent. Then, for each of our five stars, we generate a well 
populated number of test stars distributed within the errors 
shown in Figure 13] We assume a flat distribution in log T, 
with an equal number of stars on both sides from the central 
temperature, and a Gaussian distribution in logL, keeping 
the correlation between the two quantities (the diagonal er- 
ror bars) as shown in Figure p] For each of these simulated 
stars we assign masses and ages by interpolation on the grid 
of models, as well as the exact value of relative speed of evo- 
lution. Finally, we derive the best mass and age, as well as 
the uncertainties of these quantities, using weighted statis- 
tics, where the weight is the inverse of the evolutionary ra- 
pidity. The results are illustrated in Figure EJ each panel of 
which represents the mass-age distribution of each of our 
target stars color-coded according to the weights. 

We summarize our results in columns 5 and 6 of Table 
HI where we give the most probable masses and ages, and 
their associated uncertainties for our stars. The influence of 



weighting the derived stellar parameters on the evolution 
time scales is evident for stars B and E. For star B, the 
method tends to favor a slightly younger age than the central 
value of ~ 0.4 Myr shown in Figure [3] while its mass is not 
significantly affected. For star E, the large uncertainty on the 
H-R diagram leads to a correlation in the mass-age plane, 
predicting either a younger more massive, or an older less 
massive star. Nevertheless, the latter case is less probable 
than the former, leading to an estimate of the most likely 
mass of M = 13 Mq. This value is larger than that one 
would obtain by neglecting the errors in the H-R diagram, 
i.e. ~ 12 M (see Figure^. 



6 DISCUSSION ON THE DERIVED AGES 

6.1 Stellar multiplicity 

The stellar parameters derived in the previous section might 
be biased by the presence of unresolved companions. In fact, 
most of the early type stars are expected to be formed in 



binary systems (Lada 20061, with a binary fraction of the 
order of 60 - 100% ( |Zinnecker fc York^|2007| ). Given the 
distance of the Magellanic ClouddM binary systems generally 
cannot be resolved. 

As mentioned in Section HI for star B we have evidence 
from the spectrum of an unresolved later-type companion. 
For the other sources, we are unable to confirm stellar mul- 
tiplicity based on our spectroscopy. Considering that our 
spectroscopic observations for stars A, D and E were per- 
formed in two epochs, with a time baseline of ~ 2 months, 
we attempted to detect evidence of binarity through any 
radial velocity variations for these stars. We computed the 
stellar radial velocities, using the cross-correlation technique 
of |Baranne et a l. (1996). Each stellar spectrum was cross- 
correlated with a numerical template that was built from 
the stellar spectrum itself by selecting certain spectral lines. 
For this propose we rejected blended absorption lines. We 
used six absorption lines for star A and obtained an ac- 
curacy of ~2 km/s. A similar accuracy was also obtained 
for stars D and E, where we selected four absorption lines. 
While we observed no significant radial-velocity variation 
in stars A and E during the two epochs (November 2008 
and January 2009), we detected a linear trend in the radial- 
velocity variation of star D during this time window. The 
RV variation between the second (January 2009) and the 
first epoch (November 2008) is ~9 km/s. This RV variation 
can be caused by stellar activity or the presence of an unseen 
stellar companion. 

We study now in more detail how unresolved multiple 
systems can affect our derived masses and ages, and distin- 
guish two cases. 

1) The optical luminosity of the companion is signifi- 
cantly smaller than that of the primary. In this case, both 
the measured spectral type and the optical magnitude are 
representative of the primary, and thus also the position 
of the star in the H-R diagram. The presence of the com- 
panion, therefore, does not alter the derived parameters of 
the primary, but affects only the completeness of the stellar 
sample. 



1 At the distance of the LMC 1 arcsec corresponds to 50,000 AU 
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Figure 4. Illustration of the MC method to derive the stellar mass and age as well as the associated uncertainties, accounting for the 
evolution speed in the H-R diagram as a weighting quantity. For every star, the dots represent the transformation in the mass-age plane 
of the simulated population of stars within the error bars in the H-R diagram. The color is proportional (from dark to light blue) to the 
rapidity of evolution in the H-R diagram. The crosses represent the derived mean parameters, weighted according to the inverse of the 
rapidity of evolution. 



2) The optical luminosity of the companion is compara- 
ble to that of the primary. In this case, we further distinguish 
two cases, a) The two stars have identical spectral type: this 
implies that they are identical companions. They would be 
located in the H-R diagram at the same T E g we have de- 
termined through spectroscopy, and 0.3 dex fainter in logL. 
b) The two stars have different spectral types. The spectral 
type we have determined from spectroscopy will be roughly 
the average of the individual spectral types of the two stars. 
This scenario, however, would be unrealistic. This is clar- 
ified in Figure [5] where we simulate this effect on star C. 
The two companions with identical optical luminosity but 
T e ff respectively hotter and colder than what we assumed 
for star C (plotted with a star symbol) would be located in 
symmetric positions in the H-R diagram (plotted with open 
symbols) with respect to the case of identical companions 
(plotted with a filled circle) . This would lead to a significant 
difference in age between the two components, even for a 
small difference in their spectral types. Given the youth of 
LH 95, and the fact that its age is comparable to its crossing 
time (Paper I), dynamical evolution should not have altered 
significantly the binarity properties of the systems at the 
time they were formed. Moreover, since LH 95 is a loose 
system, crowding should not cause a significant source con- 
fusion. Therefore any multiple systems in LH 95 is expected 



to have be nearly coeval components. As a consequence case 
26 is not a realistic case for our stars. 

We conclude, thus, that the only plausible scenarios are 
those of cases 1 and 2a, i.e., one unseen faint component 
to our stars, or binaries with companions of nearly identical 
mass and age, each of them located in the H-R diagram 
at a luminosity 0.3 dex lower than derived in the previous 
section. Even in the scenario of case 2o, we stress that our 
derived ages would not be significantly altered, since the T e g 
of the companions would be very similar to that we have 
derived from spectroscopy. 



6.1.1 The composite nature of star B 

The issue of unresolved multiplicity is, on the other hand, 
more critical for star B, the most massive and also the 
youngest source in our sample. From our spectroscopy, we 
found clear evidence of an unresolved colder companion, be- 
cause of the presence in the spectrum of the He Ilines. The 
relative brightness of the secondary component cannot be 
negligible compared to the primary, because if this would 
be the case, the "secondary" lines would be diluted into the 
primary continuum, and not evident as they are. Since the 
only major anomaly in the spectrum is the He Ilines, we 
argue that the secondary is an early B-type star. By com- 
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Figure 5. Example of how the presence of a case 2b unresolved 
binary system, i.e., two companions of similar optical luminosity 
but different spectral types, would bias the derived position of star 
C in the H-R diagram. Instead of a single star (denoted by the 
star symbol), the two companions would be located on symmet- 
ric locations (with spectral types 07/08, 06.5/07.5, and 06/09 
respectively) at both sides from the average spectral type for this 
source (denoted by the circle symbol). The diagonal alignment 
is caused by the increasing BC with temperature. Schacrcr ct al. 
(1993) isochrones for ages 2-6 Myr are overlayed to demonstrate 
the unrealistic age differences between the companions for a bi- 
nary system (see text in Section|6.1|for details). 



paring the relative strength of the He Ilines with respect to 
the continuum, and considering that the other lines from 
the primary do not appear significantly dimmed, we esti- 
mate that the secondary (either one or more stars) accounts 
for up to 1/3 of the observed total flux in the B-band wave- 
length range. Given the very similar intrinsic (B — V) color 
for all early type stars, a 50% overestimation of the B-band 
flux results in a 50% overestimation of the V-band flux, and 
thus of the bolometric luminosity we have derived in Section 
[5] Correcting for this factor (0.17 dex in logL), the mass of 
this source decreases to ~ 60 Mq and the age reaches the 
lowest limit covered by the evolutionary models, which pre- 
dict a ZAMS massive star. 

We stress that the secondary component for star B may 
be a line-of-sight object, not necessarily physically related 
to the 02 star. Besides the colder secondary we have dis- 
cussed, we underline that even the main object might be 
a multiple system consisting of similar components. Given 
the rarity of such stars, studies of binarity for early O-type 
stars are limited to a few Galactic objects. For example, the 
02If prototype, HD 93129A is known to harbor 2 very sim- 
ilar stars (HD 93129Aa and HD 93129 Ab), re s olved thank 
to HST/FGS observations (|Nelan et aL]|2004| |20H)|). Also 



Pismis 24-1 (03.5If + 04111(f)) has been found to consist of 
at least 3 similar stars, one of the spectroscopic components 
being an eclipsing binary ( Mafz Apellaniz et al.||2007 l. For 
our star B, however, it would be hard to assume a multi- 
ple system, since the individual components would turn out 
to be even less luminous than the ZAMS (see Figure p}]), 
therefore showing parameters inconsistent with the stellar 
models. 




age (My 



Figure 6. The age distribution we derived in Paper II for LH 95 
from the analysis of the low-mass population (gaussian distribu- 
tion). The arrows, with error bars, indicate the ages and uncer- 
tainties of the five massive stars. 



6.2 Evidence for an Age Spread 

Our findings show that there is a spread in the derived ages 
for the most massive stars. This result is in agreement with 
our previous study on the age distribution of low-nrass PMS 
stars in LH 95, based on a statistical analysis of the field 
subtracted optical CMD (Paper II). There we determined 
an average age of the system of about 4 Myr, with a con- 
firmed age spread compatible with a gaussian distribution 
with a = 1.2 Myr. In Figure [5] we present the derived age 
distribution for the low-mass population, compared to the 
derived ages of the five massive members analyzed here. It 
should be noted that in this plot we preliminarily consider 
the ages obtained by neglecting any undetected binaries. 
From Figure [6] it appears that three of our sources (stars 
A, C and D) have an age which is compatible with the low- 
mass age distribution. Actually, the average age for these 
stars is exactly 4 Myr. On the other hand, star B appears to 
be somewhat younger, and star E is significantly away from 
the average population. 

Star E, however, is the faintest source of our sample, 
and thus has the largest uncertainties; moreover it is located 
in a particularly crowded sub-cluster, as shown in the stel- 
lar surface density map of Figure f7\ For this source stellar 
confusion may have contributed significantly in biasing the 
derived parameters of the star. We assess how the presence of 
unresolved equal-mass binaries, which would introduce the 
highest bias in the H-R diagram position of the stars, affects 
the observed age differences. We find that star E becomes 
even older, suggesting derived parameters which are even 
more peculiar. On the other hand, star A would turn out 
to be 1.4 Myr older, star D 1.7 Myr younger, and the age 
of star C would not be significantly affected by multiplicity. 
For these three sources, then, the average age remains com- 
patible with that of the low-mass stars even in the case of 
unresolved binarity, whereas the measured age spread would 
increase. 

Source B, as previously discussed, is the most interest- 
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Figure 7. Projected positions of our targets, in relation to the 
spatial density distribution of low-mass PMS stars in the region, 
analyzed in Paper I. The thick line delimits the central region of 
the cluster, where the density of PMS stars is above 3 standard 
deviations from the sparse background. The contours identify re- 
gions of higher concentration of PMS stars. In this study we iden- 
tified three distinct sub-clusters of such stars (see also Section 3.2 
in Paper I). 



ing case. We find that its age is younger than the others, 
at 3<r from the average age of the entire cluster, but still 
covered by the age distribution derived from the low-mass 
PMS stars. The uncertainty of this age measurement is rel- 
atively small. As we mention in Section [6. 1.1 1 this star ap- 
pears to be an unresolved composite system, and this implies 
that the stellar luminosity we have measured (see Table pf 
may be overestimated by up to 0.17 dex. This also leads 
to an even younger age than that we have reported in Sec- 
tion [5] whereas the mass of the source remains quite large 
(~ 60-70 M Q ). In Paper I we found that LH 95 is a multiple 
stellar system, consisting of three distinct sub-clusters, iden- 
tified in terms of the surface density of the low-mass PMS 
population in the region (see Figure [7b. Moreover, in Pa- 
per II we found that considering the low-mass stars, there 
is no spatial variation of either average age or age spread 
throughout the entire region of the cluster, suggesting that 
the sub-clusters are probably formed coevally. As seen in 
Figure [7] the projected location of source B is very close to 
the center of the system, but in a region characterized by 
a relatively low stellar density, and not associated to any 
of the three main sub-clusters of LH 95. Taking these into 
account, as well as its young age, it seems that source B 
was formed after the rest of the entire stellar population of 
LH 95, suggesting that the case of this star is rather peculiar. 
We can only hypothesize about the origin of a late-born 
massive star in LH 95. If source B was born in its current 
position, its formation may have been triggered by the sur- 
rounding population of massive stars. However, it is difficult 
to justify how a massive clump was prevented to collapse in 
the central part of the system for ~ 3 Myr, when massive 



star formation was taking place in the surrounding region. 
Also, the low-mass population in the immediate vicinity of 
source B does not show evidence of a younger age than the 
rest of the region. Another hypothesis is that source B has 
not formed in its present position. It is known that 10-25% of 
the O-type stars are runaways, ejected from their birthplace 
by dynamical interaction (Blaauw|1961 Poveda et al.| 1967 



Gies & Bolton 1986). Considering a typical ejection velocity 
of 40-100 km/s, and given the age of source B, of the order 
of at most ~ 0.5 Myr, it could have formed about 20-50 pc 
away from the actual position, i.e., still in the vicinity of 
the rest of the system. From our spectra, however, we do 
not measure radial velocity variations between star B and 
the other members. Also, by visual inspection on the NIR 
Two Micron All Sky Survey (2MASS) images, or the mid- 
infrared images obtained with Spitzer Space Telescope, we do 
not detect the presence of bow shocks, typically produced 
by runaway stars. Nevertheless, with no proper motion es- 
timates, we are unable to measure the actual velocity of 
this source, and therefore confirm it as a runaway star. We 
also stress that our lack of knowledge about the geometry 
of LH 95 along the line of sight limits our ability to isolate 
the plausible scenarios responsible for our results concerning 
star B. 



7 THE HIGH-MASS STELLAR IMF 

In Paper I we derived the IMF of LH 95, well into the sub- 
solar regime, based on ACS V and /-equivalent photometry. 
We found that the IMF is well approximated by a two-phase 
power law, with a break point (the "knee" of the IMF) at ~ 
1 Mq. In particular, for stellar masses larger than this point, 
up to M ~ 15 Mq, we measured an IMF slope a = 3.05, 
whereas for stars of sub-solar mass the slope decreases to 
a — 2.05. In these units, a |Salpete"r| ( |1955[ ) IMF corresponds 
to a slope a = 2.35. As mentioned before, the five stars 
studied here were not included in the analysis of Paper I, 
due to saturation in the ACS frames. In this section we 
complete the stellar IMF of LH 95, by adding the results 
of our analysis on its high-mass end. 

Since the sub-solar part of the IMF will remain un- 
changed, we limit our investigation to stars found in Pa- 
per I with masses M > 2 Mq . In this mass range the stellar 
sample inclu des about 110 members. We utili ze the method 
suggested by Maiz-Apellaniz & Ubeda ( 2005 1 to derive the 



mass distribution. These authors find that assuming a con- 
stant bin size can lead to misleading results due to the corre- 
lation between the number of stars per bin (higher for lower 
masses) and the assigned weights (from the Poisson statis- 
tics). Instead, they demonstrate that assuming a variable 
bin size, chosen such to have an equal number of sources in 
each bin, reduces the bias, which remains low also for the 
limiting case of placing a single star in each bin. Considering 
the very low number statistics in the high-mass regime in our 
sample, we choose this counting option for the construction 
of the IMF. 

We proceed by first sorting all the stellar members ac- 
cording to their mass. For every star we define the limits 
of the bin which contains only this star to be the mean (in 
logarithmic units) between the mass of the star and the pre- 
ceding and following one respectively. For the most and least 
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Figure 8. IMF for LH 95, limited to the mass interval M > 2 M Q . 
The histogram has a variable bin size, containing one star in each 
bin. The dashed line represents the best fit of the data. 



massive stars in the sample, the upper and lower boundaries 
of the bins are chosen to be symmetric from the actual mass 
of the star. The value of the mass distribution in each bin is 
then given by the number of stars per bin (1 by definition) 
divided by the width of the bin. The result is illustrated in 
Figure [8] An additional advantage of this counting method 
is that, when expressing the IMF in logarithmic units, the 
weights are identical for each source; thus we can neglect the 
errors in the fitting process of the IMF, and use a standard 
linear regression. We find that our IMF is well fitted by a 
single power low with index a = 2.32 ± 0.14, i.e., identical 
to a Salpeter IMF. The newly estimated IMF slope turns 
out to be substantially shallower than that we previously 
derived without considering the most massive stars of the 
system, included in our study here. 

In Section |6.1| we discussed how the presence of unre- 
solved stellar multiplicity may bias the derived ages for the 
massive members of LH 95. Here, similarly, we investigate if 
our measured IMF would be altered by binary systems. We 
consider the limiting case of equal mass binaries, substitut- 
ing each source with 2 stars of half luminosity, and assign 
masses by interpolation on the evolutionary models as in 
Section 5.1 For all the other sources with HST/ACS pho- 
tometry from Paper I, we assume a binary fraction / = 0.5, 
equal mass systems, and assign individual masses assum- 



ing the Siess et al. (20001 and Girardi et al. (2002) evolu- 



tionary tracks. Finally the IMF is computed and fitted as 
described above. We find under these assumptions an IMF 
slope a — 2.0 ± 0.1, therefore slightly shallower than what 
we obtained neglecting unresolved companions. This result, 
compatible with similar of such simulations in the literature 
(e.g., Sagar fc Richtler||1991 1, shows that the consideration 
of unresolved binarity makes the constructed IMF somewhat 
shallower, but not enough for its high-mass end to be signif- 
icantly affected. 



7.1 The relation between the most massive star 
and the cluster mass 

The correlation between the mass of the most massive star 



in a cluster, 



and the stellar mass of an entire clus- 



ter, M cc i, has been long investigated both through empirical 
studies and analytical modeling (e.g.,JLarson 2003 Oey & 



Clarke|2005l |Weidner fc Kroupa|200"4| [Weidner et aL|2010 | 
Its importance is fundamental in several issues, as, for exam- 
ple, it is a critical constrain to the upper mass limit for the 
formation of a star, and it enables us to investigate whether 
newborn massive stars in star clusters are randomly drawn 
from the IMF. 

In Paper I we estimated a total stellar mass for LH 95 of 
M oc i ~ 2.4 ■ 10 3 Mq. Accounting for the most massive stars 
investigated here, we have M oc i ~ 2.5 • 10 3 Mq. Consider- 
ing also unresolved binaries and a small fraction of diffuse 
low-mass stars outside the central region, we can constrain 
Moci ~ 3.0-3.5- 10 3 Mq. We utilize the model relation m max 
vs. M ec i from Weidner & Kroupa (20041, derived assum- 
ing a maximum stellar mass of 150 Mq and the [Kr oupa 



l;2001) IMF. We stress that the assumed stellar upper mass 



3-10 15 M , and m max ~ 80 M for 

O . These values are in very good agreement 

72 Mq for star B, suggesting 



limit does not affect significantly our analysis, given that 
the LH 95 cluster is not populated enough to host stellar 
masses near this limit. Moreover, as we discussed in Paper I, 
the low mass population of the region follows the Kroupa 
IMF once the unresolved binaries are accounted for, and the 
high-mass IMF derived in our study here also agrees with 
it. The | Weidner fc Kroupa| ( |2004[ ) model relation predicts a 
m max ~ 70 M© for M ecl 
M oc i = 3.5-10 3 M 
with our mass estimate of 
that the high-mass population of LH 95 is compatible with 
a random sampling of stellar masses from its IMF. 

On the other hand, our findings are in contrast with 
the empirical determinations of m max in clusters of differ- 
ent mass. As |Weidner et aL"| ( |2010[ ) point out, carrying out 
a collective research of cluster masses and maximum stellar 
mass, there is evidence that the observed m max is statisti- 
cally lower than that predicted across the entire range of 
Mod , suggesting the presence of additional mechanisms pre- 
venting the formation of the most massive stars in a cluster. 
These authors show that the relation m max vs. M oc i presents 
a plateau for clusters with M ec i in the range 10 3 - 4- 10 3 Mq, 
where the observed m max ~ 25 Mq is well below the pre- 
dicted value. Therefore LH 95 follows the empirical galactic 
relation only if star B is not a member of the system. The 
peculiarly younger age for star B seems to support this idea, 
and we have already hypothesized that this star could be the 
product of a subsequent star formation event or even a run- 



away star (see Section 6.2 1. However, both scenarios suggest 
that star B was born within the extend of the cluster, and 
therefore there is little argument against the membership of 
star B to LH 95. 



Finally, we mention that Lamb et al. (20101, by study- 



ing the O-type population of sparse clusters of the SMC, 
found examples of peculiar m max - M ec i relations. In partic- 
ular, they reported the presence of massive stars in clusters 
with a very small Meci, in strong disagreement with both 
the results of |Weidner et al.] ([2010) for the galactic clusters, 
and the predictions from the IMF. It is clear that our under- 
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Figure 9. Ha (top) and H/3 (bottom) profiles of the spectrum of 
star A. The black line represents the observed, continuum nor- 
malized spectrum. The narrow emission peaks (grey dashed line) 
at the center of the lines is the contribution from the diffuse HII 
region. The blue lines represent the expected photosphere, and 
the green lines represent the approximated line excess, after the 
photosphcric contribution is being subtracted. The red lines rep- 
resent the sum of both blue and green lines, demonstrating a good 
fit to the observed profiles. 



standing of the formation of massive stars is still somewhat 
incomplete. 



8 A B-TYPE EMISSION STAR IN LH 95 

As discussed in Section HI one star in our sample, star 
A, shows a prominent Balmer emission. Hydrogen emission 
lines in B-type stars indicate that they are either PMS stars 
still forming and emitting through accreting circumstellar 



disks (classified as Herbig Be stars; see Waters & Waelkens 



1998| for a review), or evolved B stars with hydrogen emis- 
sion, classified as classical Be stars ( |Porter fc Riviniusj2 003 ) . 
Classical Be stars are rapidly rotating B-type stars showing 
a Balmer emission produced in gaseous circumstellar 'de- 
cretion' disks ( |Porter fc Rivinius||2003"l ). Based on its T e g 
and Lboi alone this source could indeed be still in its PMS 
phase, approaching the MS. However, we exclude this pos- 
sibility for two reasons. First, the extinction of this star of 
Av = 0.88 mag (Section [5J) is very low, implying that it is 
not embedded, and second its near- and mid-infrared col- 
ors, as we present them later, are much bluer than those of 
typical Herbig Be stars ( Eiroa et al.|2002 l. 

The Hq and H/3 profiles from our combined FEROS 
spectra of star A are shown in Figure[9] The Ha line presents 
a very broad, nearly symmetric, emission that reaches rel- 
ative velocities exceeding 400 Km s" 1 . At its peak, the net 
excess is brighter than the local photospheric continuum. 
The intensity of the H/3 line is weaker in relation to the 
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Figure 10. Upper panel: the measured SED of star A based 
on our optical photometry, 2MASS, and Spitzcr/IRAC measure- 
ments (dots). The black lines represents the expected photosphere 
(T e ff = 29, 000 K) to which we applied the measured extinction 
of Ay = 0.88 mag. Lower panel: The relative ratio in the fluxes 
(observed divided by photosphere). 



continuum. Its width is comparable to that of the Ha line, 
but its profile shows two peaks. In order to characterize bet- 
ter the excess emission of star A we subtract the expected 
photospheric contribution from its spectrum. In Figure[9]we 
show also the approximate shape of the photospheric con- 
tinuum of the star (blue lines), assuming a TLUSTY ( |Lanz 
fc Hubeny||2007|) model spectrum with T cff = 29,000 K, 



logg = 4.0 and log Z/Zq = —0.5. Evidently, these two hy- 
drogen lines are expected to be in absorption from the pho- 
tosphere. We subtract the photospheric contribution from 
the observed spectrum and fit the remaining excess with an 
arbitrary function. Since the emission shows broad wings 
we find a good match assuming a line profile as the sum of 
two Gaussian distributions with equal means. The resulting 
spectrum is also shown for each hydrogen line in Figure [9] 
in green color. While our excess model is possibly oversim- 
plified, we can conclude that the observed double peaked 
shape of H/J is not real and can be well reproduced by a 
broad emission minus a narrower photospheric absorption. 
Taking the derived emission profiles into account, and given 
the absence of additional absorption from the circumstellar 
environment, we classify star A as a classical "single line" 
Be star, rather than a Shell Be star, i.e., a Be star whose 
disk is seen edge-on (Rivinius et al. 2006). 

Fast stellar rotation is a typical property of Be stars and 
it is a significant contributor to the generation of the circum- 
stellar medium surrounding these objects. The typical rota- 
tion rates of Be stars may be of several hundreds of km s _1 , 
a significant fraction of the critical rate. Our spectrum of 
star A, however, does not show significant line broadening 
due to stellar rotation. We have quantified the (projected) 
rotational velocity by considering the synthetic TLUSTY 
spectrum assumed before, convolved with line spread func- 
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tions of different values of v sin i. We have then compared 
the width of the He Iwith the observed ones. Assuming for 
simplicity rigid rotation, we obtained vsini ~ 60 km s _1 . 
This is a small value for this class of stars, and suggests that 
the stars' rotation axis is not very inclined from the line of 
sight. This is compatible with our previous argument that 
the circumstellar disk is not seen edge-on given the lack of 
Balmer self-absorption. 

Classical Be stars are also known to show an infrared ex- 
cess, commonly interpreted as free-free and free-bound emis- 
sion of the disk (|Gehrz et al.|1974| |Porter fc Rivinius|2003 



Touhami et al. 2011 1. In order to verify this characteristic for 



star A, we collected near- and mid-infrared photometry for 
this source. Specifically, we retrieved from the NASA/IPAC 
Infrared Science Archiwj(IRSA) JHK near-IR photometry 
for star A from the Two Micron All Sky Survey (2MASS) 
and mid-IR photometry in 3.6, 4.5, 5.8 and 8.0 fim from the 
Infrared Array Camera (IRAC; Fazio et al.||2004 i on-board 
the Spitzer Space Telescope ( jWerner et al.|2004| |. The result- 
ing spectral energy distribution (SED) is shown in Figure|lO} 
together with the expected photosphere after adding an ex- 
tinction of Ay — 0.88 mag, as estimated in Section [5] The 
comparison between the SED and the model photosphere 
confirms the presence of infrared excess for this source. We 
stress that the near- and mid-IR part of this SED could not 
be reproduced by any B-type photospheric model even if we 
would assume a significantly different stellar T e ff than that 
of star A. The reason lies on the fact that at these wave- 
lengths the spectrum is in the Rayleigh- Jeans regime and 
so the SED slope does not depend on T c g. Moreover, given 
the negligible effect of dust extinction at these wavelengths 
- and the fact that we do measure a very small Av - an 
inaccurate reddening estimation can neither justify the ob- 
served SED. As a consequence, the observed IR excess of 
star A is real. 

Be stars may also present line variability because of 
changes in their disk state ( jMcSwain et al.]|2009[ ). Since 
our optical spectroscopy for this source was obtained in 
two epochs (see Table [TJ, we searched for variations in the 
Balmer emission by inspecting the individual exposures from 
each epoch. Within the ~ 2 month baseline between the 
observations, we do not detect any change in the shapes, 
widths, or intensity of neither the Ha nor the H/3 emission, 
as they are shown in Figure [9] Finally, it is interesting to 
note that the inferred age of ~ 4 Myr for this star (see Table 
[2] and Figure [3j relies on the assumption that this source is 
a classical Be star, evolving off the main sequence, and not 
a forming PMS star. 



9 CONCLUSIONS 

This paper complements our investigation of the stellar OB 
association LH 95 in the LMC. In our previous studies (Pa- 
pers I and II) , based on deep HST/ACS imaging, we inves- 
tigated the intermediate- and low-mass PMS population of 
the system down to 0.4 Mq. We derived the correspond- 
ing IMF, which is found to be compatible with theKroupa 
(2001) IMF (Paper I), we determined the average cluster 
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age to be ~ 4 Myr and confirmed an age spread of the 
order of 2-4 Myr (see also Paper II). We did not find ev- 
idence of spatial variations neither in the IMF nor the age 
distribution. 

In the present study we focus on the brightest sources 
in the region, which were saturated in our previous HST 
imaging. We use ground-based high resolution optical spec- 
troscopy complemented by previous optical photometry to 
derive the stellar parameters of the five most massive stars, 
named after the letters A to E, completing thus our cen- 
sus of the LH 95 stellar population across the entire mass 
spectrum. We perform spectral classification of the sources 
in terms of the identification of specific spectral lines. We 
assign masses and ages to these stars, and we discuss the 
related uncertainties, based on comparison of their H-R di- 
agram positions with evolutionary models. We identify an 
age distribution among the stars, and we use their mass 
measurements to complete the stellar IMF of LH 95 toward 
the high-mass regime. We summarize our findings to the 
following: 

- The most massive star in the region, star B, is an 02 
giant with a mass of ~ 60 — 70 Mq. The case of this star 
is very interesting. Specifically, it appears to be somewhat 
younger (< 1 Myr) than the rest of the system. In ad- 
dition, assuming that the established relation between the 
maximum stellar mass in a cluster and the total mass of 
the cluster represents the vast majority of clusters, star B 
is more massive than expected, given the total mass of the 
stellar system (~ 3- 10 3 Mq). As such, LH 95 appears to be 
an outlier in this relation. Based on our analysis, we exclude 
any major bias that might have lead us to derive inaccurate 
stellar parameters for this source. 

- The brightest star in the sample, star A, is the second 
most massive. It is a ~ 33 Mq B0. 2 giant with significant 
hydrogen emission. Specifically, it shows very broad Ha and 
H/3 lines in emission, originating from a gaseous circumstel- 
lar disk, as well as clear near- and mid-infrared excess. Con- 
sidering its low extinction and blue infrared colors, we assess 
that star A is not a Herbig Be PMS star evolving toward the 
MS, but a classical Be star, evolving away from it. 

- The average ages and age spread of the most massive 
stars are consistent with those previously derived by us for 
the low-mass PMS populations. Exception to this behavior 
is star E, which appears significantly older than the rest of 
the stars. We assign this peculiarity to the severe crowding 
by low-mass stars in its immediate surroundings that affects 
its photometric measurements, and its low brightness, which 
introduces large uncertainties in our measurements. 

- The complete stellar IMF of LH 95, from the most 
massive stars down to the intermediate mass regime with 
<; 2 Mq follows a Salpeter slope. Together with the low- 
mass stars, therefore, this IMF is compatible with a Kroupa 
IMF. 
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